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Abstract. We use the long-wavelength formalism to compute the bispectral non-Gaussianity 
produced in two-field inflation. We find an exact result that is used as the basis of numerical 
studies, and an explicit analytical slow-roll expression for several classes of potentials that 
gives insight into the origin and importance of the various contributions to /nl- We also 
discuss the momentum dependence of /nl- Based on these results we find a simple model that 
produces a relatively large non-Gaussianity. We show that the long-wavelength formalism 
is a viable alternative to the standard 5N formalism, and can be preferable to it in certain 
situations. 
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1 Introduction 



During inflation the energy density of the Universe is assumed to be dominated by the 
potential energy of one or more scalar fields in order to have a sufficiently rapid expansion to 
solve the homogeneity problems that plagued pre-infiationary cosmology (horizon, flatness, 
etc.). Very importantly, inflation also provides the initial adiabatic density perturbation that 
generated the large-scale structure observed today. Observations of the fluctuations in the 
cosmic microwave background radiation (CMB), in particular those made by the WMAP 
satellite, have verified the basic predictions of infiation. The most important observational 
parameters so far from the point of view of infiation have been the amplitude and the slope 
(spectral index) of the primordial power spectrum, as well as some limits on the amount of 
tensor perturbations and the running of the spectral index. 

Unfortunately this small number of observational parameters means that a very large 
number of quite different infiation models are all still consistent with the data. To further 
narrow down the number of viable inflation models additional observables are required. A 
very promising candidate is the non-Gaussianity of the primordial power spectrum. In its 
most simple form this is encoded as a non-zero three-point correlator of the CMB temper- 
ature fluctuations, or equivalently a non-zero bispectrum, which is the Fourier (or spherical 
harmonic on the sphere) transform of the three-point correlator. The quantity defined as the 
bispectrum divided by the power spectrum squared is called /nl- The current limits on this 
parameter /nl (assuming momentum dependence of the local type, relevant for this paper) 
after seven years of WMAP data are —10 < /nl < 74 at 95% confidence level [1]. The newly 
launched Planck satellite [2] is expected to significantly improve these constraints, down to 
la error bars of about 3-5 (depending on the use of polarization data) [3, 4]. While standard 
single-field slow-roll inflation predicts an unobservably small value of /nl [5], many other 
models predict much larger values that could be detected or ruled out by Planck. 

Both supersymmetric particle theory and string theory suggest the existence of multiple 
scalar fields that can infiuence the early Universe. If more than a single scalar field plays a 
role during inflation, isocurvature fluctuations will be produced in addition to the adiabatic 
one. While these isocurvature fluctuations might have directly observable consequences in 
the CMB [1] , in this paper we are only interested in the effect of the isocurvature fluctuations 
on the adiabatic one during inflation. This effect can be important even if the isocurvature 
fluctuations disappear after inflation. The important point here is that while in single-fleld 
inflation the adiabatic perturbation is constant on super-horizon scales, this is no longer 
true in multiple-fleld inflation. In fact, the isocurvature perturbation acts as a source for 
the adiabatic perturbation on super-horizon scales and this source is multiplied by the i]-^ 
parameter [6]. This r/-*- (defined properly in the next section) is proportional to the component 
of the field acceleration perpendicular to the fleld velocity. In other words, rj-^ is non-zero if 
the field trajectory makes a turn in field space. Only during such a turn will the isocurvature 
mode influence the adiabatic one on super- horizon scales (see also [7]).^ 

There are two main ways to produce non-Gaussianity during inflation: during and after 
horizon crossing of a perturbation mode. Horizon crossing is defined as the moment when 
the physical wavelength (k/a)^^ of the fluctuation becomes equal to the Hubble (or horizon) 
length H~^, i.e. when k = aH (with k the wave number of the mode, a the scale factor 

^This last statement is strictly only true on a flat field manifold (trivial field metric) with standard kinetic 
terms. On a curved manifold 77^ can be non-zero even for a straight field trajectory because of the connection 
terms in the covariant derivatives, see [8] and [9]. In this paper a trivial field metric will be assumed. 
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of the Universe and H = a/a the Hubble parameter). The first type of non-Gaussianity 
is produced in all inflation models, but it is unobservably small (i.e. slow-roll suppressed), 
unless the model contains non-standard kinetic terms (higher derivatives), like DBI inflation 
[10, 11] for example. In this paper we will not consider those models, but instead focus on 
the super-horizon type of non-Gaussianity. As is clear from the previous paragraph, super- 
horizon non-Gaussianity can only be produced in multiple-field inflation models where the 
field trajectory makes a turn in field space. To compute this type of non-Gaussianity we will 
make use of the long-wavelength formalism developed by Rigopoulos, Shellard and Van Tent 
[6, 12, 13], hereafter refered to as RSvT. 

The main purpose of this paper is twofold. In the first place we want to further work 
out, simplify and study the general analytic expression for /nl of RSvT in the case of two 
fields only. In the second place we want to clarify a few remaining formal issues with the 
formalism, and compare with an alternative formalism for computing /nl, called the 6N 
formalism [14-18]. 

As will be shown, the expression for /nl simplifies significantly, although a final inte- 
gral remains that cannot be done analytically. As a starting point for numerical work this 
expression is very useful though, and it gives a fully exact numerical result since no slow-roll 
approximation is used for the super-horizon evolution. However, to be able to derive explicit 
analytic results, we have to make use of the slow-roll approximation. Using this approxima- 
tion we derive explicit analytic results for a number of classes of infiationary potentials. Some 
of these have been treated before in the literature using the 5N formalism, but others are 
new. In particular we show that models with a potential of the form W{(j), a) = acjf + f3a'^ 
with p = q will never give a large /nl that persists until the end of infiation (unless infiation 
somehow ends right during the turn of the field trajectory, but in that case a very careful 
treatment of the transition at the end of inflation will be required) . We also present a simple 
model that does produce a "large" /nl of the order of a few. The reason we choose this 
model is that it can be treated not only numerically, but also analytically. 

The formalism of RSvT also allows us to compute the momentum dependence of /nl 
due to the fact that different modes cross the horizon at different times. This effect has 
usually been ignored in the literature where it was often assumed that /nl is momentum- 
independent (see e.g. [19]), although recently people have started looking into this [20, 21]. 
Here we compute this momentum dependence in an exact way and show that, depending on 
the model, it can lead to relative effects of order 10% even within the range of momenta that 
are observable by Planck. 

We have extended the formalism of RSvT with an exact treatment of the second-order 
source term at horizon crossing. While negligibly small in the models we consider, the 
inclusion of this term allows for an exact analytic comparison with the 6N formalism. We find 
that our analytic slow-roll results agree exactly with those derived using the 5N formalism, 
where available. For models where slow roll breaks down long after horizon crossing and 
which have to be treated numerically we also find excellent agreement. 

Apart from providing an alternative way of computing the bispectrum, which is always 
useful, the long-wavelength formalism provides a number of advantages compared to the 6N 
formalism. Very importantly, the long-wavelength formalism allows for a simple physical 
interpretation of the different terms, showing the contributions from adiabatic and isocurva- 
ture modes and making clear why some of them can become big and others cannot. While 
we do not pursue this in the present paper, the formalism also provides the solution for 
the second-order isocurvature perturbation and hence the isocurvature bispectrum could be 
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computed as easily as the adiabatic one. 

Many people have worked on non-Gaussianity, both predictions from inflation and es- 
timators for CMB observations. The reduced bispectrum has been given in [22] for the 
equilateral type, in [3, 4] for the local type and in [23] for the orthogonal type of the pri- 
mordial bispectrum. The different shapes were studied in detail in [24, 25]. Bispectrum 
estimators were developed in [22, 26-29]. All kinds of inflationary models have been stud- 
ied as well. For instance, we have learned that single-field inflation models cannot produce 
large non-Gaussianity [5], unless some non-trivial potential is used [30] or higher derivative 
contributions are introduced as for the Dirac-Born-Infeld action [10, 31-33] or K-inflation 
[34, 35]. The study of the effective theory of inflation has also turned out to be very fruitful 
[36]. These models were extended to incorporate multiple fields in [11, 32, 37-39]. Large 
non-Gaussianity can also be produced at the end of inflation [40-45] or after inflation, in 
models with varying inflaton decay rate [46] and in curvaton models [47-52]. 

Large scale evolution of perturbations during inflation up to second order became pos- 
sible through their consistent gauge invariant definition [53-55]. Within the SN formalism 
several authors have investigated the bispectra of specific multiple field inflation models [56- 
61]. Two-field models, being easier to deal with, have gained popularity though. Vernizzi 
and Wands studied the double field sum potential [19], while the double product potential 
was studied in [62]. Conditions for large non-Gaussianity were found in [63]. 

The paper is organized as follows. In section 2 we summarize the long-wavelength 
formalism of RSvT and define the various quantities used in the paper. Here we also describe 
the second-order source term, which is a new extension of the formalism. At the end of the 
section we give a brief overview of the 6N formalism for comparison purposes. In section 3 
we work out the general expression for /nl in the case of two fields. We also compute the 
momentum dependence that arises in the case that not all scales cross the horizon at the same 
time. These expressions, derived without using any super-horizon slow-roll approximation, 
are one of the main results of the paper and the starting point for our numerical analyses. In 
order to find completely explicit analytic expressions for /nLj however, we do need to assume 
slow roll, as well as some conditions on the potential. This is treated in section 4, where we 
also compare our analytic results with those obtained using the 6N formalism, for as far as 
the latter exist. In section 5 we use the two-field quadratic potential in order to compare our 
exact numerical results with those of the (5A^-formalism. We also present a simple potential 
that can produce an /nl of the order of a few, which falls into the category of potentials that 
can be treated analytically, thus allowing us to test our results. We conclude in section 6. 
Finally, in the appendices we give supplementary information on the basis in field space 
that we use, compute the second-order source term, comment on some gauge and formal 
issues, and provide several intermediate steps of our calculations. Note that A introduces in 
particular a small improvement of the basis defined in [8] that makes it more convenient for 
numerical calculations during periods when the fields oscillate. 

2 Basic equations and definitions 

This section sets up the starting point for the work in the following sections. It is mostly 
a summary of the long- wavelength formalism and its results as presented in [6], although 
the part about the second-order source term in section 2.1 and some results on the Green's 
functions in section 2.2 are new. Section 2.1 describes the non-linear equations for the per- 
turbations, section 2.2 shows how to solve them using Green's functions, and section 2.3 gives 
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the formal expressions for the two and three point correlation functions of the perturbations, 
or rather their Fourier transforms, the power spectrum and the bispectrum. It is the latter 
that will be worked out in great detail in the rest of the paper. Finally in section 2.4 we give 
a brief overview of the 6N formalism for comparison purposes later in the paper. 

2.1 Non-linear equations 

In the long-wavelength formalism space-time is described by the long-wavelength metric [64] 



where the lapse function N{t, x) defines the time slicing and a{t, x) is the space dependent 
scale factor. The Hubble parameter is defined as H = dtlna/N. In order to simplify 
superhorizon calculations we choose to work in a flat gauge NH = 1, that is we choose time 
slices in which the expansion of the universe is homogeneous and the time variable coincides 
with the number of e-folds t = In a. 

On the matter side we assume in this paper two scalar fields with a trivial field metric, 
although the formalism can in principle deal with an arbitrary number of scalar fields living on 
an arbitrary field manifold. The energy-momentum tensor for the two fields (f>'^ {A, B = 1,2) 
is 



where W is the potential. The Einstein summation convention is assumed throughout this 
paper. We also define the derivative of the fields with respect to proper time as = dt(t)^/N, 
with length 11. 

The long-wavelength formalism corresponds to the leading-order approximation of the 
spatial gradient expansion (see [65, 66] and references therein). In this expansion all quantities 
are expanded in terms of a small parameter 1/{HL), where L is the characteristic physical 
length scale of the perturbations (i.e. proportional to a). The leading-order approximation 
of the spatial gradient expansion is equivalent to neglecting the k'^ term (which comes from 
the second-order spatial gradient and is of order 0{1/{HL)'^)) with respect to the 0{\) terms 
in the equation for the perturbation modes. Because of the very rapid growth of a during 
inflation, this is in principle a well-justified approximation from just a few e-folds after horizon 
crossing of the perturbation mode under consideration, when the decaying mode will have 
disappeared. However, as pointed out in [67, 68], if slow roll is broken at horizon crossing and 
for some e-folds afterwards, a cancellation of the 0{1) terms can cause the decaying mode 
to remain important during this period. In those papers it was shown that for single-field 
inflation there may be an enhancement of the curvature perturbation both at first and second 
order due to the effect of the fc^ term even on super-horizon scales, if the decaying mode has 
not yet vanished. In this paper we will assume slow roll to hold around horizon crossing, 
so that the decaying mode will quickly disappear, and the long-wavelength approximation is 
valid on super- horizon scales. 

The Einstein and field equations are given by [12] 



ds^ = -N'^{t, x)dt^ + a^{t, x)dx'^ 



(2.1) 




(2.2) 




H = 



2H 



-SH^- 



diH 



,2 



(2.3) 



H 
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with = SttG = STr/rripi and W^a = dW/dcf)^. While the first three equations look hke 
background equations, they are actually fully non- linear since H{t,x), (l)^{t,x), and U^{t,x) 
defined above are functions of both time and space. We define an orthonormal basis in 
field space (see A for details) through the field velocity, = n^^/n, and successively higher- 
order time derivatives of the fields [8]. This basis allows us to easily distinguish effectively 
single-field effects with m = 1 from truly multiple-field effects with m >2. The local slow-roll 
parameters then take the form 

e{t,x) = 7]^t,x) = ri-^{t,x) - ^"^^^ 



H' ' ' n ' ' ^ ' ' n ' 

Xit, x)^W22 + e + , it, x) ^ - er?ll , {t, x) = - er?^ , (2.4) 

where Wmn = Wm,n/(3-ff^) and Wmn = G-^G^W^AB- Throughout this paper the indices 
l,m,n will indicate components in the basis defined above, taking the values 1 and 2. The 
correction term in the expressions for .^H and comes about because the proper definition 
is = {[{I / N)dt]'^Il'^) / (H'^U) and NH = 1. We note here that we have not made any slow- 
roll approximations so far; the above quantities should be viewed as short-hand notation and 
can be large. We also give the time derivatives of the slow-roll parameters, 

e = 2e(e + r/ll), j)" = + (r/^)^ + (e - t)^ = + (e - 2r?H)?7^, 

X = er?!! + 2ex - (r/")' + 3(r/^)2 + + ^r?^^^ + VF221, (2.5) 

and of the unit vectors, 

= v^el 4 = -V^et. (2.6) 



Here we have defined Wimn = {V^e/ K)Wimn/{'^H'^)- Hence the tilde has a different meaning 
in the case of two and of three derivatives, but it is these specific combinations that always 
appear in the equations. 

Our main variable describing the perturbations is [12] 

C = Smidi In a - ^ [en^AdiC^^] , (2.7) 

with m = 1 the adiabatic component and m = 2 the isocurvature one. The non-linear 
quantity has been constructed to transform as a scalar under changes of time slicing on 
long wavelengths. In the single-field case and when linearized it is just the spatial gradient of 
the well-known curvature perturbation Of course in the gauge we have chosen 5^ Ina = 
and the first term in the expression for (^^^ disappears. The exact evolution equations for 
and its time derivative 6^ = dt {CT) [6] 



Vj, 



,{t, x) + Aab{t, x)vib{t, x) = 0, where Via = (C^ C^^ of) , {21 



so a, 6 = 1, 2, 3. We have simplified the system by omitting the time derivative of the adiabatic 
component, dj, since it is given by 9} = 2r]-^(f, valid fully non- linearly [6]. The matrix Aab 
in the two-field case is 

■ -2r/-L 

-1 I . (2.9) 

3x + 2e2 4er/l! + 4(r/-L)2 + ^11 3 + e + 2r/ll 
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To solve this equation, we expand the system as an infinite hierarchy of hnear inhomo- 
geneous perturbation equations. To first and second order we obtain 

v'i^ + A^(t)v'^ = li^{t,x), (2.10) 
+ ^1? W-!? = -A^\t, x)vil^ + bfjit, X). (2.11) 

Here the source terms 6^^-* and ft*-^-* have been added to describe the influence of the short- 
wavelength modes on the long-wavelength system given in (2.8), providing the necessary 

is found by perturbing the exact A matrix, giving A^^^ 



initial conditions. A^^^ is found by perturbing the exact A matrix, giving A^^^{t,x) = 
A^^) (t)v^\t,x). The explicit form of A is given in (D.6), where we have dropped the 



superscript (0) for notational convenience. We have also defined v^^ = 

The source term b^^^ can be expressed in terms of the linear mode function solutions 

x2rl, using a window function yV{k) which guarantees that short wavelengths are cut out to 
get only the contribution to the long- wavelength system, 

hf^ = I j0j^Mk)Xillik)al{k)ik,^^^^^ + c.c. (2.12) 

The quantum creation (aln) and conjugate annihilation (a^) operators satisfy the usual 
commutation relations. The linear mode solutions can be determined exactly numerically, or 
analytically within the slow-roll approximation (which, as observations indicate, seems to be 
a very good approximation at horizon crossing). See section 2.2 for explicit expressions for 
these x2rn as well as for the window function W. 

The source term 6)^ was either neglected by RSvT [6] because it is small, or in earlier 
papers (e.g. [13]) approximated by perturbing xim, which turned out not to be a good 
approximation. While this contribution to /nl is indeed small, here we compute it explicitly 
in order to allow for an exact comparison with known results in the literature. We find that 
it can be expressed by means of the window function as 

/ (|^/(|^mmax(.',.)) 

x|L,,,(^)x(i;(fe^^)x«(fc,^)a3„(fc04(fc)i(fc^ + fc^)e'('^'+')■^ 

+^^a^,c(^)^^i](A:^^)X«(fc,^)a1„(/c')at(fc)ifc,e^('^'+^)•^ + c.c.|, (2.13) 
where the derivative of the window function peaks at the scale that exits the horizon last. 

(2) 

We have split b^J into a local part proportional to Labc and a non-local part proportional to 
Nabc- In order to find these factors, we generalize the results of Maldacena [5] to multiple 
fields. This is similar to the work done in [56], but due to the different definitions and gauge 
choices used in that paper, we found it easier to rederive the results from scratch. Maldacena 
computed the third-order action for in the uniform energy density gauge. In order to 
calculate the three-point correlation function he performed a redefinition of C to remove 
terms in the action which are proportional to the equations of motion. Our generalization of 
this calculation can be found in B, while the explicit expressions for the components of Labc 
and Nabc are given in section 2.2. 
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2.2 Green's functions 

Equations (2.10) and (2.11), together with the initial condition Via {t — t- — oo) = 0, can be 
solved using a simple Green's function Gab{t,t')- In matrix notation it satisfies [6, 13] 

^G(t,t') + A(t)G(t,t') = 0, G{t,t) = l. (2.14) 

Starting from this equation, to lighten the notation, when we write A we actually mean 
i.e. the matrix in (2.9) with all local slow-roll parameters replaced by their background version 
that depends on time only. Looking at this equation of motion and its initial condition, we 
see that the solution can be written as 

G{t,t') = F{t)F-\t'), (2.15) 

where F(t) satisfies the same equation of motion (2.14) as G{t,t') with an arbitrary initial 
condition. From this we immediately derive that 

^G{t,t')-G{t,t')A{t') = 0. (2.16) 

The solution of (2.10) and (2.11) can now be written as the time integral of Gab con- 
tracted with the terms on the right-hand side of these equations: 



v'il\t,x) = J -A^^^^(A,,t)a1„(fc)iA;,ei^- + C.C., (2.17) 



with 



Vam{k,t)= / dt'Gab{t,t')W{k,t')xl2{k,t') (2.18) 

J —QO 

and 

vll\t,x) = - f dt'Gab{t,t')Abcd{t')v^^Xt'^x>d\t'^^) + f dt'Gab(.t,t')b'i\t',x). (2.19) 

J ~oo J—oo 

As before, v^^^ = d~'^d^vfj . 

Written in components (2.14) gives the following equations for the Green's functions: 

^Gi,.(t,t') = 2r?^(t)G2.(t,0, 

^G2..(t,0 = G3x(t,0, (2.20) 
^G3..(t,t') = -A;2{t)G2.{t,t') - A;-i{t)G^^{t,t'), 

Gab{t,t) = 6ab- 

We can also rewrite this as a second-order difi'erential equation for G2x' 

^G2At, t') + A33(t)^G2.(t, t') + ^32(t)G2.(t, t') = 0. (2.21) 
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For the derivatives with respect to t' we find: 



^G^2{t,t') = -2r]^{t')6^i+A32{t')G.3{t,t'), 

^G.3(t,t') = -G.2{t,t')+A33{t')G^3{t,t'). (2.22) 

The solutions for the x = 1 components of (2.20) are simple: Gn = 1, G21 = G31 = 0. 
To find the solutions for the x = 2, 3 components we assume that we have found a solution 
g{t) that satisfies (2.21). Then a second, independent, solution is given by 

dfyffl, Yit) . ^e-/ = ^^J^y (2.23) 



Hence 



and 

r u /^ ^(*) r a ^ 3it)Yit) g{t) g{t') g{t)Y{t) 

Of course Gr3{t,t') = 2 dt 77^(^)^23 (i, t') and Gi2(t,t') = 2 //, dt r/^(t)G22(t, t')- For exact 
calculations the Green's functions will be determined numerically, but in an approximate 
slow-roll treatment we can sometimes find analytic solutions, see section 4. 

For the linear mode solutions at horizon crossing, X^, we will assume in this paper 
the analytic slow-roll solutions determined in [8] . Observations of the spectral index indicate 
that slow roll is a good approximation at horizon crossing. Note however that, with the 
exception of section 4, we do not assume slow roll to hold after horizon crossing. Moreover, 
the assumption of slow roll at horizon crossing is not a requirement to compute these linear 
solutions, we could just as well numerically compute the linear mode solutions exactly. For 
the window function used in the calculation of the linear solution we take a step function, 
see [6, 12], so that its time derivative is a delta function: W = 6{kc/ {aH\/2) — 1), where c is 
a constant of the order of a few, e.g. c = 3. Then 

Vamit) = Gabit,t,)X^^{t,)Qit-t,), (2.27) 

where the step function Q{x) equals 1 for x > and for x < 0. The time t=K is defined 
by aH = kc/V2, i.e. a time slightly after horizon crossing when we have entered the long- 
wavelength regime. While results right at of course depend on the details of the window 
function, a few e-folds later any dependence on W has disappeared. Moreover, under the 
assumption of slow roll at horizon crossing, all quantities change very little between horizon 
crossing and t*, so that final results do not depend on the choice of c and can be taken 
equal to the horizon-crossing time determined from k = aH in the final expressions. Defining 
7* as 7* = —KH^/{2k^/'^.ye^), where the subscript * means evaluation at t = t^,, the matrix 
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X.^^\U) is given by x[Wu) = X2^{t^:) = 7*, X^'^{t^:) = —x*l*, the other components being 
zero [6]. Hence we have 

vn = 7*G(t - t^), vi2{t) = 7^, {Guit, U) - x*Gi3{t, U)) Q{t - U), 

V21 = 0, V22it) = 7, (G22(t, U) - X*G23{t, Q) Q{t - Q, 

V31=0, V32{t)=-f*iG32{t,U)-X*G33{t,U))Q{t-U). (2.28) 

We also define the short-hand notation Vam by Vam{t) = 7*0(t — t^:)vam{t)- 

For the second-order horizon-crossing solutions we find from (B.IO) and (B.ll) (see B) 
that the slow-roll matrices Labc and Nabc have elements satisfying 

-^111* = £*+??!!, -^^122* = — ( e* + f]* ~ Xi 



I* ; 



-^^211* — V* ) -^^222* — V* 

Lll2* + -^^121 = "^V*, ^112* + Ni2U = -"^llii 

L212* + -^>22i* = 2 (e* + r/l - x*j , N212* + ^''221* = X*, (2.29) 



with the other elements of Nabc being zero. As explained in the appendix, a slow-roll ap- 
proximation which expresses 6'^ in terms of has been used. This means in particular that 
the subscripts a, b, c only take the values 1 and 2, but not 3. However, for consistency in the 
notation, we will define here all entries of L^bc and N^bc to be equal to zero if one or more of 
the indices are equal to 3. 



2.3 Two and three point statistics 

So far we have used time slices on which the expansion of the universe is homogeneous 
(5jlna = 0), since it simplifies super- horizon calculations. However, to make contact with 
the proper gauge-invariant expression for C it turns out to be necessary to change to uniform 
energy density time slices {dip = 0). On such slices the adiabatic perturbation variable has 
the simple form 

Cl = dilna. (2.30) 

At first and second order the relation between the adiabatic component of Q in the two 
gauges is (see [6] and the discussion in C) 

C^^t) = Cf (2.31) 
Cp(0 = Cp(t) + 2r?^C(^)^Cp, (2.32) 

where again 

^(1)1 = d-'^dxl^^\ Indeed one can show that not only do we end up with a 
total gradient through this gauge transformation, we also obtain the gauge-invariant quantity 
corresponding to the curvature perturbation ^: in the flat gauge and for superhorizon scales 
one can show that [69] (see also [55] for the energy density definition of Q) 

^(2)1 = 5^^(2)1 _ ^(1)1^(1)1. (2.33) 

The second term on the right-hand side cancels exactly the gauge transformation term and 
we are left with the space gradient of the gauge-invariant quantity. 
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From (2.17) we see that the Fourier coefficients of C^^^^{x,t) are given by Cj^^^i't) = 
vim{k,t){aln{k) + am{—k)). Hence the power spectrum, that is the two-point correlator of 
the Fourier coefficients, is 

{CS^& = '^'('^i + k2)vim{ki,t)vim{h,t). (2.34) 

Conventionally (see e.g. [70, 71]) a quantity V(^, also called the power spectrum of (, is 
defined to remove the overall delta function and the factor coming from the vim (see 
(2.28)), as follows: 

^3 

^C(^' ^) = ^ i)vim{k, t). (2.35) 

The scalar spectral index is then deffiied as 

„ _ 1 ^ dl^n ^ dlnPg dt. ^ dlnV^ 1 

d In A; dt* dink dt* 1 — e* ' 

where we used that k = aH\/2/c and H = —eH for our gauge where t = lna. 

By combining the different permutations of {C^i^ C^^^ Ck^^) Fourier components 

of the linear and second-order adiabatic solutions (first subtracting the average of Q'^^^{x, t) 
to get rid of the divergent part), we find the bispectrum^ [6] 

iCkClCl/''' = {27:)-^lH\Y.k,)[f{krM) + f{kiM) + f{k2M)] 

s 

= {2TT)-^/H\Y.K)B(^{kiMM), (2.37) 

s 

where 

/(/C, A;') = Vim{k)vin{k') ^ ri-^V2m{k)vin{k') + ^Gia{t, tk')Mabc{tk')Xbm{k, ty)Xcn{k' , tk') 

-^J dt'Gia{t,t')AabcVbm{k)Vcn{k')j +k^k', (2.38) 

where k' refers to the scale that exits the horizon last and Mate = Late + ^^abc- Finally 
we introduce the parameter /nl, basically defined as the bispectrum divided by the power 
spectrum squared, which gives a relative measure of the importance of non-Gaussianities of 
the bispectral type (see e.g. [5, 19]): 

6. _ Bc^{ki,k2,k^) 
:/NL 



/(fci,A:2) + /(fcl,fc3) + /(fc2,fc3) 

{ki)vin{k2)vin{k2) + 2 perms. 



%Vi^{ki)%Vi^{k2) + {k2 o A;3) + {ki o A;3) 

(2.39) 



The quotient is called — |/nl and not simply /nl because it was originally defined in terms 
of the gravitational potential $ and not C as $ = $l + /nl ('^'l ~ (*^l)) t^]- During recom- 
bination (matter domination) the two are related by C, = — |*^*- Moreover, when computing 



^In the literature (e.g. [5]) one often sees a factor in front of the bispectrum (as well as in front of 

the power spectrum). This is due to a different definition of the Fourier transform. We use the convention 



where both the Fourier transform and its inverse have a factor (27r) 



-3/2 
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the bispectrum divided by the three permutations of the power spectrum squared using this 
expression of ^ one obtains 2/nl due to the two ways the two $l inside the second-order so- 
lution can be combined with the two hnear solutions to create the power spectrum. Together 
these two effects explain the factor —6/5.'^ 



2.4 JA^- formalism 

An alternative formalism to compute /nl is the so-called (5A^-formalism [14-18] . In order to 
compare our results of the next sections to those obtained using the 6N formalism, for those 
cases where the latter are available, we give here a brief overview. 

The 6N formalism uses the fact that the adiabatic perturbation (^^ on large scales is equal 
to the perturbation of the number of e-folds 5N{t,tif) between an initial flat hypersurface at 
t = t^,, which is usually taken to be the horizon crossing time, and a final uniform density 
hypersurface at t. One can then expand the number of e-folds in terms of the perturbations 
of the fields and their momenta on the initial flat hypersurface 

Fin f)N ^ FfiN 

So instead of integrating the evolution of through equations (2.10) and (2.11) one can 
evaluate the derivatives of the number of e-folds at horizon crossing and thus calculate C^. 

Because of the computational difficulty associated with the derivatives with respect to 
n^^, slow roll is assumed at horizon exit so that the terms involving the momentum of the 
fields can be ignored. This is a crucial assumption for the 6N formalism. The final formula 
then reads 

f)N ^ rfiN 

mt,t.)^^Ht + -,^^HtH?. (2.41) 

up to second order. From it one finds the following expression for the bispectrum: 

iCkClCl)^'^ = lNANBNcD{6K6cl,g{6<p'' * 5</.^)fe3) + perms., (2.42) 

where * denotes a convolution and the average of {6(f) * S4>) has been subtracted to avoid 
divergences. N^a denotes the derivative of N with respect to the field 4>'^. Using Wick's 
theorem this can be rewritten as products of two-point correlation functions to yield finally 

— — /NL <5N = ^"^5~- (2.43) 

Notice that this result is momentum independent and local in real space, although attempts 
to generalize to a scale-dependent situation have recently been made in [20, 21]. This formula 
can be used numerically or analytically to calculate /nl- However, for any analytical results 
and insight one must assume the slow-roll approximation to hold at all times after horizon 
exit (see for example [19, 59, 62]), except for the special case of a separable Hubble parameter 
[72, 73]. 



^In [6] and earlier papers we used a slightly different definition of /nl which was larger by a factor of -18/5. 
Here we conform to the definition that is now generally accepted in the literature. 
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3 General analytic expression for /nl for two fields 



In this section we will further work out the exact long-wavelength expression for /nl for the 
case of two fields, given in (2.39). No slow-roll approximation is used on super- horizon scales 
in this section. In particular this means the formalism can deal with sharp turns in the field 
trajectory after horizon crossing during which slow roll temporarily breaks down. In the first 
subsection we restrict ourselves to the case where ki = k2 = to lighten the notation. In the 
second subsection we show how the result for /nl changes in the case of arbitrary momenta. 

3.1 Equal momenta 

In the case of equal momenta, equation (2.39) reduces to 



6^^^ _ -Vlmit)vinit) 
it)vim{t)) 



/nl = . ,;.2 S / dt'Gia{t,t')Aabc{t'Hm{t'Wn{t') " 27]^ {t)v2m{t)vinit) 



oo 



-Gla{t,U)Mabc*Vbm.{U)Vcn{t*) )■ (3.1) 



We remind the reader that indices l,m,n take the values 1 and 2 (components in the two- 
field basis), while indices a,b,c,... take the values 1, 2, and 3 (labeling the C^, C^j and 
9'^ components). To make the expressions a bit shorter, we will drop the time arguments 
inside the integrals, but remember that for the Green's functions the integration variable 
is the second argument. Using the result (D.l) proved in D.l we can write Aabi as a time 
derivative and do an integration by parts, with the result 



dt'GlaAabcVbmVcn = 2r]-^V2mVln+ / dt' Aab-T-j[GlaVbmVln] + / dt'GlaAabcVbmVcn, (3.2) 

50 J —OO J —OO 



am 



-oo J —oo J —oo 

where the index c does not take the value 1. Here we used that the linear solutions v, 
are zero at t = —oo (by definition), that the Green's function Gia{t,t) = 6ia, and that 
Alb = —2ri^6b2 (exact). We see that the first term on the right-hand side exactly cancels 
with the gauge correction (the second term in (3.1)) that is necessary to create a properly 
gauge-invariant second-order result. 

We start by working out the second term on the right-hand side of (3.2). We find 

/* d /"* d 

= Aab*Gia{t, U)vbm*Vln* + 7*0(* " **) dt' Aab [GldAdaVbmVln - GiaAbdVdmVln 

— GiaVbmAldVdn] 

= Aab*Giait, t^)vbm*Vln* " 7*0(* " **) dt' AabAidGiaVbmVdn, (3.3) 

where, as before, a subscript * means that a quantity is evaluated at f*. Using the explicit 
form of the matrix A (2.9) and the solutions Vam (2.28) this becomes 



/ = %@{t - t*)6m2Snl [-2rit + X*Gl2{t,t^) + A32*Gl3{t,U) - X*A33*Gl3{t,U)j (3.4) 
+jl@it - t^)6jn2Sn2 dt'2r}^V22 -2-r]^ V22 - G12V32 + A32G13V22 + A33G13VZ2 ■ 
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From now on we will drop the overall step function, which just encodes the obvious condition 
that t > t^,. Realizing that A32V22 + ^33^32 = —^^^32 we can do an integration by parts: 



I = ll^m2 



5ni ( -27]^ + x*Gi2{t, U) + ^32*Gi3(t, U) - x*Mz*Gi^{t, U] 



-Sn2'^V*X*Gl3it,Q 



(3.5) 



+7* ^m2S: 



'n2 



2V {V22r + ( -2Gi2 + A33G13 + ^Gl3 ) ^'22?^32 + ^13(^32) 



To this result we have to add the final term on the right-hand side of (3.2). Using the explicit 
expression for the matrix A and doing some more integrations by parts this can be worked 
out further, as can be found in D.2. The final result for /nl in the equal momenta limit is 
(including also the final term of (3.1)) 



6 



/nl 



12 



, /_ xono 1 9iso + 9sr + 9int 

[1 + {vi2yY \ 



(3.6) 



where 



(e + i^){v22f + ^^22^-32, 



e^+ril r]^vi2 3 / || 



9int 



dt' 



2(?? ) (f22) + (e + ??")^22^^32 + (^'32) - Gi3l'22 Cw22 + 9/? V32 



■ (3.7) 



Here we have defined 



C = 12r]^x - Qvh^ + 6(r/ll)2r/^ + 6{r]^f - 2r/^e" " 2r?l'e^ - ^(VFan + W^222), (3.^ 



where as before Wimn = {V2e/ K)Wimn/{3H'^). We should add that although no slow-roll 
approximation has been used on super-horizon scales, we did assume slow roll to hold at 
horizon crossing, in order to use the analytic linear short-wavelength solutions (2.28) and to 
remove any dependence on the window function W. Observations of the scalar spectral index 
seem to indicate that slow roll is a good approximation at horizon crossing. In a numerical 
treatment we could use the exact numerical solutions instead. 

Looking at (3.6), which is one of the main results of this paper, we can draw a number 
of important conclusions. In the first place there is a part of /nl, namely the first term in 
Qsr, that survives in the single- field limit. It corresponds to the single- field non-Gaussianity 

(2) 

produced at horizon crossing and comes from the source term. It agrees with the single- 
field result of Maldacena [5] for f^^. The rest of the result is proportional to ui2, which 
describes the contribution of the isocurvature mode to the adiabatic mode. In the single- 
field case it is identically zero, so that there is no super-horizon contribution to /nl in that 
case. Moreover, since 6^ = 2rj^C,'^ , such a contribution only builds up when rj^ is non-zero, 
i.e. when the field trajectory makes a turn in field space. We also see that there are three 
diff'erent sorts of terms in the expression for /nl • The gsr terms are proportional to a slow-roll 
parameter evaluated at and thus are always small because we assume slow roll to hold at 



- 14 - 



horizon-crossing. Altliougli the terms proportional to V12 and l/fi2 in gsr are time varying, 
one can easily show that neither ^12/(1 + ^^12)^ nor 7)^2/ (1 + ^12)^ ^^^^ bigger than 0.33. 
The giso terms are proportional to V22, the pure isocurvature mode. These terms can be 
big, in particular during a turn in field space, but in the models that we consider, where 
the isocurvature mode has disappeared by the end of inflation, they become zero again and 
cannot lead to observable non-Gaussianities. The reason that we do not consider models with 
surviving isocurvature modes is that in that case the evolution after inflation is not clear. In 
the presence of isocurvature modes the adiabatic mode is not necessarily constant (indeed, 
that is the source of the non-Gaussianities we are considering here), which means that the 
final results at recombination might depend on the details of the transition at the end of 
inflation and of (p)reheating. Hence we will make sure that in all models we consider the 
isocurvature modes have disappeared by the end of inflation, which means in particular that 
the turn of the trajectory in field space has to occur a sufficient number of e-folds before the 
end of inflation. Note, however, that this is a constraint we impose voluntarily to simplify 
the evolution after inflation, it is in no way a necessary condition for the validity of our 
formalism during inflation. Finally, the third type of term in (3.6) is the integral in gint- It 
is from this integrated effect that any large, persistent non-Gaussianity originates. 

For completeness we also calculate the power spectrum, which according to equa- 
tion (2.35) takes the simple form 

and the spectral index, calculated analytically using equations (2.36), (2.16) and (2.5), 



— 1 



1 — e* 



4e, - 2r]l + 2-^^ ( - 2rit + X*vi2 



+Gi3{t,U) (-^"221* + 2e^ + r/f + 77,^2 + 3e*(ryU - X*) - 2??Ux* + 



(3.10) 



3.2 General momenta 

We turn now to the more general case where each scale exits the horizon at a different time 
tfc. , defined by aH = kic/V2, where c ~ 3 is a constant allowing for some time to pass 
after horizon exit so that the long-wavelength approximation is valid (see the discussion 
in section 2.2). It is important to realize that it is not the momentum dependence of the 
bispectrum that we are discussing here, but of /nl- The momentum dependence of the local 
bispectrum is dominated by the momentum dependence of the power spectrum squared, 
leading to the well-known result (see e.g. [24]) that it peaks on squeezed triangles where one 
of the momenta is much smaller than the other two. Here we are discussing the momentum 
dependence of /nl, so one has divided by the power spectrum squared. This /nl, often called 
/^L in the 6N literature, is usually assumed to be momentum- independent. However, as we 
will show this is not true and its momentum dependence can lead to relative effects of order 
10% even within the range of momenta that are observable by Planck. 
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Assuming fci > ^2, i-e. t^j > ^^2) find that (2.38) reduces to 



f{ki,k2 



7li7fc2 



dt'Giait, t') [AahAie + AaeAn,] ^{,^^2 (^O^enfci (*') 

+Aab{tki)Gla{t, tk 1 ) ['^bmk2 {tki)hn + Vimk2{'tki)Sbn] 
ki ) y^cmk2 ^bmk2 



(3.11) 



where again we have used the result (D.l) and have done an integration by parts that 
cancels the gauge correction term as in (3.2). The indices h and h do not take the value 
1. We have introduced the notation vnki = 5ii and Vi2ki = Gj2(t, ifcj — XkiGi3{t,tki), where 
Xki is evaluated at t^, . We notice that due to the step functions the integral's lower limit 
corresponds to the time when both scales have entered the long-wavelength regime, i.e. the 
time when the larger ki (smaller wavelength) exits the horizon. The expression has become 
more complicated as compared to (3.2) and (3.3) since the Vbmki refer to a different initial 
value depending on the horizon crossing time of each scale ki . Following the same procedure 
as in the previous section we find that 



where 



6„ , , . f{kl,k2) + f{k2,k3) + f{kuk3) 

^ 7fci74[l + ivi2k,y][l + {vi2k2r] + 2 perms. 



(3.12) 



fiki,k2) = -27^^ 7^2(^12) giso{ki,k2) +gsr{ki,k2) +gint{ki,k2) +gkiki,k2) (3.13) 



with 

giso{kl,k2) = (e + r/ll)(v22)^ +W22t'32, 
G22kik2Vl2ki 1 



gsr{kl,k2)=Vki 



G22kik2 N \ , ^Xk2 

H : — L^33A:ifc2 



2 fi2fe2 2wi2fci / 4 

3, ,XkJ^vx2k^ ^ \ ^k,+r]l^ 

"2 '•^^^i + VkJ^22kik2 + \ + <^22fcifc2 



gint{kl, k2) 

gk{ki, k2) 



dt' 

1 



4wi2fci 



Wl2fc2 ^ " J 2(i;i2)2 

2(r/-^)^(t;22)^ + (e + r]^^)v22V32 + (v32f -Gi3V22{Gv22 + 9V'^V32) 

'iGi3{xkiG22kik2-Xk2G33kik2) + G32kik2 + +27y| - Vi2ki^ 



+ ^Gi2A:ifc2(-2??4 +XkiVi2ki) + ^G32kik2{VkiGl3 " 1) 
Gi2kik2 ( , II I ± - \ 1 / ,11 Xkx \ ^ 

for ki > k2 > k3 and C was defined in (3.8). We introduced the notation 

(Vu)"^ = Vi2kiVi2k2, {V22f' = V22k^V22k2, {^32)'^ = V32k^V32k2, 

?^22^^32 = \{v22kiV32k2 + ?^22fc2 ^32A:i ) , 



(3.14) 
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and also Gjjfcifc2 = G'jj (t^^ , tfc2 ) > while the subscript on the slow-roll parameters denotes 
evaluation at the relevant time that the scale exits the horizon. The Green's functions that 
appear without arguments denote G{t,tk^) outside or G{t,t') inside the integral. 

Although this expression is quite a bit longer than (3.6), there are many similarities 
between the two results. The whole expression is again proportional to Vi2ki, except for the 
single-field horizon-crossing result, so that there is no super- horizon contribution to /nl for 
the single-field case. In the giso and gsr terms we recognize the familiar terms of the equal- 
momenta case, i.e. the isocurvature contributions proportional to V22ki as well as the horizon 
crossing terms now evaluated at t^^ and (note that for ki = k2, Gnk^ki = 1 identically 
and we regain the expressions of (3.7)). The integral has also retained its form. The rest of 
the terms, namely those in g^, are terms arising due to the different horizon-crossing times 
of the scales and are identically zero for the equal-momenta case ki = k2 where Gijk-^k2 — ^ij- 
All terms inside gk are proportional to a slow-roll parameter evaluated at horizon crossing 
(using the fact that G13 = 6*12/3 up to slow-roll corrections, see (4.2)), except for the very 
last term on the second line. However, 032^1^2 is expected to be quite small: for ki = k2 
it is zero, and for ki ^ k2 it becomes the linear solution for the isocurvature velocity 9'^ 
(see (2.28)). Hence we do not expect gk to give a large contribution, which is confirmed 
numerically. As we will see later on, the dominant contribution to the differences between 
different momentum configurations comes from the changes in the other terms. 



4 Slow-roll approximation 

While the exact result for /nl, equation (3.6) or (3.12), is an extremely useful starting point 
for an exact numerical treatment, the integral cannot be done analytically. In order to find 
explicit analytic results that will be very useful to gain insight and draw generic conclusions, 
we need to simplify the problem by making the slow-roll approximation. In subsection 4.1 
we further work out (3.6) under this approximation. Even then the integral can only be 
done analytically for certain specific classes of infiationary potentials, which are treated in 
the other subsections. 

4.1 General expressions 

Considering the slow-roll version of equation (2.21) we find that g{t) (as defined above equa- 
tion (2.23)) satisfies 



We see that Y{t) oc exp(— 3t) so that f{t) is a rapidly decaying solution that can be neglected 
(see (2.23) for definitions). After the decaying mode has vanished the solutions for the Green's 
functions simplify to 



9 + X9 = 0. 



(4.1) 



G32{t,t') 




9{t') ' 

-X{t)G22{t,t'), 




(4.2) 
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4.1.1 Equal momenta 

Using the last two relations in (4.2) and dropping higher-order terms in slow roll, (3.7) reduces 
to 



giso = (e + r?ll -x){v22f 

lb 



_ + -ql r]fvi2 _ 3 
2vl^ + 2 2 



e* + r/H - X* + ^ ) , (4.3) 



9int 



{V22 



+ G12 ( v^x-'2vh^--{W2ii + W222] 



Inserting these terms into (3.6) we find an expression that can be considered the final expres- 
sion for /nl in the slow-roll approximation, and is the one that will be used in section 4.3. 
It also proves useful, however, to rewrite it in a different way using integration by parts. 

We use the slow-roU version of relation (2.22), 2r]-^ = —-^Gi2{t,t') + xG\2{t,t'), to do 
an integration by parts, leading to 

{e^ + r]l - x*)x* 



gint = Vi2 + 



+ / dt'Gi2{v22y' 



(e-Fr?ll -x)x^ 
2r/-L 



-27711 ry^ - ^(1^211 + W222) + + 



± , (e + 7?ll -x)x 
2r/-L 



(4.4) 



Using the slow- roll version of the relations (D.ll), 

='ieri\\ +{rl\f + {r]^f -Will and = 3er/^ + 2r/llr/^ - r/^X " ^211, (4.5) 

as well as the time derivatives of the slow-roll parameters in (2.5), we can derive that 



dt 



e + rjW - x)x 
2?7-L 



1 

2rp 



+ (e + j?ll)x^-4feryil +(r?^) 



X 



(4.6) 



+4 {e^ri\\ + e{rj\f - ei-q^f + ri\\{r]^f^-{e + rj^ - xWiii 



2r?^ + 



(e + r?ll -x)x 



W2U + (e + - 2x)W^22i 



Inserting this into expression (4.4) for gint and including the remaining terms in the expression 
for /nl we finally obtain 



6 -2(^12)^ 

5 [1 + {vuyr 



(e + ^ll_;^)(^22)' ' '^*^'' ■ (^*+^*-^*)^* 



I ^ — ^ 
2vl, + 2 



V12 



I II I ^* 

V12 



+ 



f dt'Gl2{v22) 



-X+e+T/" 



(r7-^)2\ 1 ~ - Y - 

+o(^211-VF222-4^22l) 



^ + V" -X 
2r/-L 



X 



1^111-^221 r 1^211 



(4.7) 
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This is the alternative final result for /nl in the slow-roll approximation. 

Equation (4.7), as well as (4.3), is characterized by the same features as the result of the 
exact formalism. We can easily distinguish the pure isocurvature V22 term, which we assume 
to vanish before the end of inflation in order for the adiabatic mode to be constant after 
inflation, as well as the terms evaluated at the time of horizon crossing, which are expected 
to be small. Any remaining non-Gaussianity at recombination has to originate from the 
integral. In subsections 4.2 and 4.3 we will further work out the expressions of this section 
for the case of certain classes of potentials to gain insight into their non-Gaussian properties. 
But first we look at the momentum dependence of /nl hi section 4.1.2. 

4.1.2 Squeezed limit 

In this section we will calculate the slow-roll expression for /nl in the case where k = k-^ <^ 
ki = k2 = k' , what is usually refered to as the squeezed limit. Note that we assume ki = k2 
for simplicity, to keep the expressions manageable, it is not a necessary condition. We start 
with equation (3.12) and follow the procedure of the previous subsection, that is we use the 
slow-roll approximations (4.2) for the Green's functions and drop higher-order terms in slow 
roll. Since there are only two relevant scales the expression simplifies to give 



6 



-2vi2k'/[l + ivi2k'?] 

'^[l + {vi2k'y] + 2[l + {vuk) 



l'^vi2k'\giso{k', k')+gsr{k', k')+gint{k', k') 



+2vi2k ygtso 



{k' ,k)+gsr{k' ,k)+ gintik' , k)+gk{k', k) 



(4.8) 



where 7 = 7/^/ / 7^ and 
giso{k',k) = (e + ??ll - x)v22kV22k 

G22k'kVl2k' 




2 v^fc' + Vk')G22k'k+ ^ 




+ G22k' 



efc' + %' 



gintik', k) 
gk{k',k) 



dt'v22kV22k' 



Vl2k 



2r/^ -r/^+- 



2r/J 



2vi2kVl2k' ' 

+ G'i2(r/^X-2r/llr/^- 



-(^"211 + ^222) 



-^{Xk'G22k'k — XkG33k'k) + Y^G32k'k{ 

I2k'k -Vk' 7. ^ + 

\ 2 V Vl2k' 



+G 



-6 + ek>+2rfl, + 2riyVi2k' 



I "Xk' /-, , - 
+ -^{^ + Vl2k' 



(4.9) 



The first line of (4.8), proportional to 7^, comes from the f{k',k') term and it is identical 
to expression (4.3). The difference is that now it occurs with a weight 7^ compared to the 
terms originating from f(k',k) that come with a weight 2. Obviously, in the case of equal 
momenta where Gijkik2 = ^iji the expression reduces to equation (4.3). 

The 7 terms can be safely neglected in the squeezed limit because 7^ scales as e"'^'^*, 
where At is the number of e-folds between horizon exit of the two scales. If for example 
the two scales exit the horizon with a delay At ~ 7, which corresponds to k' ~ 1000/c, 
approximately the resolution of the Planck satellite, we find that 7^ ~ 10~^. 
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The functions vi2 and V22 increase and decrease respectively (from their initial values 
and 1) only a little until the turning of the fields. The later the relevant scale exits the 
horizon, the less time there is available for Vi2 to evolve, so the smaller is the value of V12 
(and the larger for V22) during this period. During the turning of the fields isocurvature 
effects turn on. Both V12 and V22 vary wildly during this period. V12 grows and reaches a 
constant value afterwards, while V22 varies and reaches zero when isocurvature effects cease. 
In the models we studied we found that while during this period V22 continues to behave in 
the same way, i.e. being larger for the scale that exits last, V12 changes behaviour and also 
becomes larger for the scale that exits last. In the end we observe that /nl in the squeezed 
limit is smaller than in the equal-momenta case that was treated in the previous subsection. 
The effect is particularly pronounced during the turn of the field trajectory, mainly due to 
Qiso- As we will show in section 5.1, these effects can reduce the value of /nl during the turn 
of the field trajectory by 10% on scales that are within the resolution of Planck. 

4.2 Potentials with equal powers 
4.2.1 Quadratic potential 

The quadratic potential has been widely examined in the past and it is known that it cannot 
produce large non-Gaussianity (see for example [19]). Here we use our results to analytically 
explain why. While the quadratic potential is a special case of the more general sum potential 
treated later on, it is still interesting to discuss it separately in a different way. We start by 
deriving the result that for a quadratic two-field potential within slow roll. 

Working out the right-hand side, using (2.5), we find 

X = 2e + 7?ll-^-\ + ^. (4.11) 

Inserting the relations (4.5) (with the third derivatives of the potential equal to zero, since 
we have a quadratic potential) this reduces to 

X = 6 + ??ll-^. (4.12) 
77II 

It can be checked that this result does indeed satisfy the general equation for the time 
derivative of x (2-5) within the approximations made, and the remaining integration constant 
is fixed by realizing that this result has the proper limit in the single- field case. This concludes 
the proof of (4.10). 

Since the third-order potential derivatives as well as the first term of the integral in (4.7) 
are identically zero, we find that for a quadratic potential the integral completely vanishes in 
the slow-roll approximation and no persistent large non-Gaussianity is produced. Numerically 
we find that even for large mass ratios, when during the turn of the field trajectory slow roll 
is broken, the integral is still approximately zero, see section 5.1. 

Using this result (4.10) for x we can also solve (4.1): 

m = ^, (4.13) 
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and hence find that 

Note that even though g{t) is a large quantity, of order inverse slow roll, it is still slowly 
varying, as we have shown, with its time derivative an order of slow roll smaller. 

4.2.2 Potentials of the form W = a(j)P + (3aP 

For a potential of the form 

W{(t),a) = a4>P + ^a" (4.15) 

we can work out explicitly the form of the integrand in equation (4.7). We have to use the 
slow-roll version of equations (2.3) and (2.4) to easily find after substitution that 

Qint = — A o -at, (4.16) 

Ju 2k'^ {a<j)P + l3<jPyf{a{p-l)(l)Pa'^-l3y{py-l)<i,'^aPyf 

where y = q/p. 

From this expression we can derive an important result: for y = 1, i.e. p = q, we 
immediately see that the integral is zero. This means that no persistent non-Gaussianity 
can be produced after horizon exit for potentials of the form W{<l), a) = acjP + /5o"P, at least 
within the slow-roll approximation. This generalizes the result for the two-field quadratic 
potential of the previous subsection to any potential with two equal powers. 

4.3 Other integrable forms of potentials 

In general, the first step of finding an analytical expression for the integral gint is to solve 
the differential equation (4.1) for g in order to determine the Green's functions. To do that, 
one tries to express x as a time derivative of some other quantity. In the slow-roll limit 

Wii = e-r,\ W2i = -7i^, (4.17) 

so that X can be written as 

X = 2e + W22-Wu. (4.18) 

which, as can be checked, cannot be expressed as a derivative of a known quantity for a 
general potential. Thus we are forced to examine special classes of potentials. 

4.3.1 Product potentials 

First we consider potentials of the form 

W{(l),a) = U{ct))V{a), (4.19) 

inspired by the analytical study done in [62, 63]. From our point of view, the advantage 
of these potentials is that their mixed second derivative Wipa can be expressed in terms of 
the first derivatives to finally give for the second-order derivatives of the potential in the 
adiabatic and isocurvature directions: 



W21 = [W4,^- W,, + 2e(ef, - e\^) ei<^ei,, (4.20) 
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where we used (A. 8) to eliminate the unit vector 62 in terms of ei. It is straightforward to 
show that the second-order derivatives in the directions of the basis vectors are related: 



2e + W22 - Wii _ eia _ei^ 
W21 eio- 



(4.21) 



so that only two of them are independent. Now we can use (2.4) and the above results to 
write X 



x = w, 



21 



eio- _ ei^ 



dt 



ln(ei0eiCT) , 



(4.22) 



where the derivatives of the unit vectors are given in (2.6). Hence looking at equation (4.1) 
we can identify the Green's function g to be 

g{t) = ei^{t)ei„{t). (4.23) 

After a few more manipulations the integrand of Gi2{t,t') in (4.2) takes the form 

ld5 



4 dt 



(4.24) 



where S = ef^ — ef^, so that the analytical form of the two independent linear perturbation 
solutions in the slow-roll approximation is 



V12 



2ci(^*ci(T> 



V22 



(4.25) 



The final step is to write the integrand of gmt in (4.3) in terms of the potential's 
derivatives and rearrange terms to form time derivatives. One can prove that then the 
integrand can be rewritten as 



1 



9int 



dt' 



S{t) - S{t')] (W^^{t')ej^{t') - W^^{t')el{t') 



dt'. 



(4.26) 



After performing the integration and adding the rest of the terms we find that 
_6 _2{S^S^nSl-ll 

g/NL — (1 + S'^ 2SS )^ ^9iso ~r gsr "r gint) j 



(4.27) 



where now 



5^-1 



9int 



2{S - S,) 

52-1 



l|\ 1 + 3S{S - 2S,) + 2s: 



e* -I- 77* 



+ ryil - x* 



2 -3 + g2 + 4SS^ - 2Sl 
^* 2K 



51+1 
252 



(4.28) 



Comparing to the results of [62, 63] we find complete agreement. 



*Note added: Very recently the same results were obtained in [74] for the transfer functions Tns and Tss 
of product and sum potentials, which turn out to coincide with V12 and V22- 
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Looking at the result for /nl for the product potential we can draw a number of con- 
clusions. The only time-dependent slow-roll parameters appear in giso- These terms and 
consequently /nl can vary significantly during a turn of the field trajectory but, as explained 
before, in the models we consider isocurvature modes have disappeared by the end of inflation 
so that the adiabatic mode will be constant after inflation, which means giso will disappear 
again and cannot give any persistent non-Gaussianity. The rest of the terms involve slow-roll 
parameters evaluated at horizon crossing, which are small. Hence we conclude that any large 
non-Gaussianity will have to come from the denominator becoming very small (since IS'I < 1 
the numerator cannot become large) to compensate for the small slow-roll parameters at 
horizon crossing. We see that this can only happen when 6*, 5* ^ ±1. In the remainder of 
this section we will study the two different cases that satisfy this condition: a 90° turn in the 
field trajectory {S = — S*), or the same field dominating both at the beginning and at the 
end (S* = 5*). 

First we study the case where the field trajectory makes a 90° turn. The field (f) is 
dominant right after horizon crossing, which means |eio-*| <C 1, |ei^*| 1 and hence 5* 
1. Later on occurs a turn in the field trajectory and afterwards a leads inflation, so that 
\^14>\ ^ 1) kicrl ~ 1 and S* — t- — 1. Then we find that both gsr and g^so go to zero, which 
means in particular that we satisfy the condition on the disappearance of the isocurvature 
mode that allows us to directly extrapolate the results at the end of inflation to the time of 
recombination. The non-zero term comes as expected from g^nt and it is given by: 

- ;g/NL = e^+r}l-x* = -W^a*, (4.29) 

since Wo-a* = W22* = X* ~ ^* ~ V*- Hence we see that for any product potential where the 
field trajectory makes a 90° turn no significant non-Gaussianity will be produced, at least 
within the slow- roll assumptions used to derive this analytic result. 

Next we look at the opposite limit, where one of the fields, (j), is dominant both at 
horizon crossing and at the end of inflation. This means |eicr*| <C 1, |ei^*| ^ 1 and |eio-| <C 1, 
|ei</,| ~ 1, so that S"* — t- 1 and 5—7-1. This includes the case where we have a perfectly 
straight field trajectory, i.e. an effectively single-field situation, where obviously no super- 
horizon non-Gaussianity is produced. However, we find that even more generally in this limit 
the contributions from g^so and gint go to zero and we are left with only the single-field result 
from gsr- 

-^/NL = e* + r?l (4.30) 

Hence no significant non-Gaussianity is produced in this limit. 

We conclude that if we impose the condition of the disappearance of the isocurvature 
mode by the end of inflation, to simplify the evolution afterwards, the product potential can 
never give large non-Gaussianity, at least within the slow-roll approximation. 

4.3.2 Potentials of the form VF((^, a) = {U{(j)) + V{a)y 
Next we consider potentials of the form 

WicP,a) = iU{<j>) + Via)r, (4.31) 
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which, to our knowledge, have not been worked out before for general u.^ While of course not 
the most general two-field potential, it can accommodate potentials with coupling terms of 
the form a^(/>^ + /3^cr^ + 2al3(j)a or higher-order combinations. Note that in the case of = 1 
the potential becomes the simple sum potential, which has been studied before [19, 63]. 

Just as for the product potential, we find that mixed second derivatives of the potential 
can be expressed in terms of the other derivatives: 

9 ~ 9 4e(i/ - l)ehej^ _ , _ Ae{iy - 1)6^ ,ef^ 



2e{v - l)ei4,eia{el„ - ej^) 



(4.32) 



Again there are only two independent second derivatives of the potential in our basis: 



21 



22 



(4.33) 



'Icr ^10 

Following the procedure of the previous section we rewrite x a-s 

ei^\ d 



2e 

X = — + W2i 



ei0 eio- 



dt 



In H^/' 



ei(/>eio- I , 



and then find an analytical expression for g, 

g{t) = F2/-(t)ei^(t)ei,(t). 
The integrand of Gi2{t,t') is now written as 



1 fK 



2\ ^1^ 



dZ 
d^' 



where Z = Ve^^ — Ue\^. Finally we find that 



V12 



Z — Zjf 



^^22 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



Rewriting the integrand of gint in terms of the potential's derivatives yields after a few 
manipulations 



9int 



w; 



-2/u 



p2 p2 



d 

At' 



• d ^ W'/''{t')e{t')e\^{t')eUt') 



Wy^it' f^'^ ^^^'\ w..{t')e%{t') - Ty^^(t')eL(iO)] |(4.38) 



and adding the rest of the terms results in 
6 



NL- 



2wr''{z-z.)\l^A.. 



2 I 9iso ~l~ Qsr ~\~ Qint 



(4.39) 



^Note added: While we were doing the final editing of our manuscript, a paper [75] appeared on the arXiv 
where the authors studied this type of potential using the 5N formalism. Their result for /nl agrees with 
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where 




i(f>* 



'la*/ 



0^2/^^ 2 2 



W^/i^ 2 2 



9int 



e 

+ - 




(4.40) 



Note that the first term on the second hne of gint is also related to the pure isocurvature 
mode, but we have not incorporated it in giso in order to remind the reader that it originates 
from the integral. 

As in the case of the product potential we will study two limiting cases, to get some 
insight into the behaviour of /nl- First is the limit where the field trajectory makes a 90° 



turn. We assume that (j) dominates inflation at horizon exit, that is |eio-*| ^ 1, |e 



10* 



1 and 



— >■ 14. At late times, after the turn of the field trajectory, the second field a is dominant 
and the contribution of (j) is negligible, so that |ei0| <C 1, |eio-| ~ 1 and Z — )• —U. Then we 
find that gsr and giso go to zero, while the remaining contribution to /nl comes from gint, as 
expected, 

6 „ [/, + K U^ + V^ f ,11 



5 



U + V, 



e^+rjH - X* 



(4.41) 



So we see that we need a significant decrease in U between horizon crossing and the end of 
inflation, as well as a relatively small value of K, to get a large /nl- Of course we cannot 
increase U-t/U too much without breaking slow roll. In section 5.2 we investigate numerically 
the properties of a model with a sum potential and confirm the validity of the above limit. 

In the opposite limit (j) dominates both at horizon crossing and at the end of inflation, 
i.e. |eio-*| <C 1, \ei^^\ ^ 1 and |eio-| <C 1, |ei^| 1 so that Z^ — )■ 14 and Z ^ V. Then the 
expression reduces to 

- ^/nl = - y_y X*) , (4.42) 

which comes from gint- Note that we have assumed here that V ^ V^:. In the (effectively) 
single-field case this is not valid; in that case we find that gint and giso are zero and gsr goes to 
the single- field result, e* We remark that in this limit giso is zero, so that the adiabatic 

mode is conserved after inflation. In order to make /nl large, one might be tempted to take 
V close to 14. However, that means a does not evolve and we are in an effectively single- field 
situation, where the above limit is not valid. Instead the situation is somewhat similar to 
the previous limit: we need a large value of [/* and relatively small values of 14 and V to 
overcome the small values of the slow-roll parameters at horizon crossing. It might not be 
simple to satisfy these conditions together with the requirements of this limiting case that (p 
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Figure 1. The time evolution of the fields (left) and the field trajectory (right) during the period of 
the turn of the field trajectory, for the model (5.1) with initial conditions 0o — ctq = 13k^^ and mass 
ratio 7710 /mcr = 20. 



dominates both at horizon crossing and at the end of inflation, with a period of a domination 
in between; we did not further study those types of models. 

As a final remark we point out that the power v of the potential does not appear 
explicitly in the limits for /nl- Of course its value will play a role in determining the field 
trajectory and the values of the slow-roll parameters, but that is only a relatively small effect. 
We have verified this result numerically for several values of the power v of sum potentials 
of the form (4.15). 

5 Numerical results 

The formalism we have developed so far provides a tool to calculate the exact amount of 
non-Gaussianity produced during inflation driven by a general two-field potential, beyond 
the slow-roll approximation. While we assumed slow roll in the previous section, in order to 
derive analytical results, we return here to the exact formalism for a numerical treatment. 
In the following subsections we investigate the properties of the quadratic potential as well 
as a potential of the sum type that can produce an /nl of the order of a few, and compare 
our results to those of the 5A^-formalism. 

5.1 Comparison with 6N for the quadratic potential 

We investigate the quadratic potential 

W = \ml4>^ + ]^mW (5.1) 

choosing our parameters as follows: m^lmfj = 20, ttt-o- = IO^^k^^ and the initial conditions 
00 = o'o = 13k^^ at t = for a total of about 85 e- folds of inflation. From now on we will 
denote the heavy field as cf). We choose to present this particular mass ratio because the fields 
oscillate wildly during the turn and slow roll is badly broken, so that it provides a serious 
check both of our formalism and the 6N one. Of course we have also run tests with smaller 
mass ratios when slow roll is unbroken and verified our analytical slow-roll results. 

We solve the field equations (2.3) numerically and in figures 1 and 2 we plot the values 
of the fields, the unit vectors, and the slow-roll parameters as a function of time during 



- 26 - 



1.0 



0.5 



0.0 



-0.5 



-1.0 





ei<r 










15 


- 








e 




— ^10 










10 
5 
















\ 

\ 










-5 














\ 
I 
1 


/• N \ / 
' ^ \ / 

' * vy 












\ 1 \' 
\l 








/ 


1 > _ 






-10 








'il 




41 


42 43 


44 


45 46 


47 


48 




41 


42 


43 


44 45 46 


47 48 



N 



N 



Figure 2. The unit vectors (left) and the slow-roU parameters e,ry" and rj-^ (right) as a function of 
time during the turn of the field trajectory, for the same model as in figure 1. 



the range of e-folds where the heavy field <j) is approaching zero and starts oscillating. In 
the beginning of inflation (j) dominates the expansion while rolling down its potential and 
about 40 e-folds after the initial time t = it starts oscillating around the minimum of its 
potential. The heavier (p is, the more persistent are the damped oscillations. During the 
period of oscillations the unit vectors, as well as the slow-roll parameters e, t?", and rj-^, 
oscillate too. For m^/rriu = 20 the maxima of the slow-roll parameters are much larger than 
unity and slow roll is temporarily broken. During these oscillations the light field a starts 
driving inflation and rolls down its potential until it also reaches its minimum and starts 
oscillating. We take the end of inflation when e = 1 during this second period of oscillations. 
The situation is similar to the limiting case we studied in section 4.3.2 with |eio-*| ^ 1 and 
\ei4>\ <. 1. 

In figure 3 we plot the /nl parameter as calculated in our formalism both the numerical 
exact version (3.6) and the slow-roll analytical approximation (4.7) (but using the exact 
background), as well as the result computed numerically in the context of the (5 A'^- formalism. 
The horizon-crossing time is defined as 60 e-folds before the end of inflation. We do not 
expect any large non-Gaussianity to be produced in this model, since we have shown that 
the integral of (4.7) is equal to zero in the slow-roll approximation. The final value of /nl 
calculated in all three cases is O(10~^). Our results coincide completely with those of the 
(5A^-formalism, thus reinforcing the validity of both formalisms. We also show /nl for a much 
smaller mass ratio, m(f,/ma = 4, where slow roll remains valid throughout the turn of the 
field trajectory, verifying our analytical slow-roll result. 

The peak of the /nl parameter during the turning of the fields is due to the isocurva- 
ture terms giso in the slow-roll analytical formula. As expected this effect is transient and 
disappears when the isocurvature mode V22 has been fully converted to the adiabatic one. 
There is no surviving isocurvature mode in this model. The higher is the mass ratio, the 
larger is the magnitude of the peak as a consequence of the more violent oscillations. 

For completeness, we plot in figure 4 the power spectrum (3.9) and the spectral index 
(3.10) of this model. We see there is a jump in both of them during the oscillatory period of 
the heavy field, but afterwards they become constant again. 

Finally in figure 5 we plot the exact numerical /nl in the squeezed limit and in the 
equal momenta limit. As mentioned in section 4.1.2, we see that the /nl parameter in the 
squeezed limit is smaller than in the equal-momenta one. From figure 5 we can see that for 
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Figure 3. We plot the /nl parameter for the model (5.1) with initial conditions (j)o = = 13k~^ 
and mass ratio m^/m^ = 20 (left) and m^/mg. — 4 (right). The red line is the exact numerical result, 
while the blue dot-dashed line shows the slow-roll analytical approximation (but using the exact 
background). We also show the numerical 6N result as the black dashed line, which lies practically 
on top of our red result. 
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Figure 4. The exact numerical power spectrum (left) and the spectral index (right) for the same 
model as in figure 1. 



k' = 1000 A; (roughly corresponding to the Planck resolution) the peak value of /nl is more 
than 10% smaller than for k' = k, for this particular model. 

5.2 A simple model producing large non-Gaussianity 

In this section we introduce a model that produces an /nl of the order of a few, which is 
two orders of magnitude larger than the single- field slow- roll result. So in that sense we can 
call it large. From the point of view of observations with the Planck satellite it is probably 
still a little bit too small, but we have taken this particular model to be able to make the 
connection with our analytical results. 

The /nl limit (4.41) that we calculated in section 4.3.2 can be simplified for the sum 
potential = 1) to give 

" ^^^^ = -kWTWV ^^-^^ 
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Figure 5. In the first plot we depict the /nl parameter for the equal momenta limit in red and the 
squeezed limit result for fc'/fc = 1000 in dashed black. In the second plot we show the dependence of 
the discrepancy of the first peak value of /nl in the two limits on the ratio k' /k. We used the same 
model as in figure 1 and t^i = 25 (60 e- folds before the end of inflation). 




Figure 6. The time evolution of the fields (left) and the field trajectory (right), for the model (5.3) 
with initial conditions (fiQ = 18K~^,fTo = O.OIk"^ and parameters a2 — 20k~^,62 = 7k^^ and 64 = 2. 
Only the time interval during the turn of the field trajectory is shown. 



where we used the definition of Wmn and the slow-roll version of (2.3) for H. We can easily 
infer that in order to obtain a large value for /nl, the heavy field (j) should end up with a small 
value at the end of inflation, while a should obey a potential characterized by a large second 
derivative and a small value at horizon crossing. Such properties can be accommodated by 
a potential of the form 

V{a) = bo- b2a^ + b^a^ (5.3) 

with Bq = 62/(4^4) so that the minimum of the potential has W = U + V = 0. 

To illustrate the above we investigate a model with 02 = 20k~^,62 = 7k^^, 64 = 2, and 
initial conditions 0o = 18k~^ and ao = 0.01k~^, so that the light field is standing on the local 
maximum of its potential, for a total amount of 85 e-folds of inflation. This type of effective 
potential might be realized in the early universe during second-order phase transitions. We 
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Figure 7. The unit vectors (left) and the slow-roU parameters e,?]" and r]-^ (right) as a function of 
time during the turn of the field trajectory, for the same model as in figure 6. 



solve the field equations (2.3) numerically and in figures 6 and 7 we plot the evolution of the 
fields and the unit vectors, as well as the slow-roll parameters. The situation is qualitatively 
the same as in the case of the quadratic potential: in the beginning (j) dominates inflation 
while rolling down its potential, then there is a period of violent oscillations around (j) = 0, 
and a takes over and starts rolling down towards the minimum of its potential. 

The behaviour of the unit vectors is that of the limiting case we studied in section 
4.3.2, that is \eia*\ ^ 1 and \ei^\ <C 1. We will try to obtain an analytical estimate for the 
magnitude of /nl- The final value of /nl is reached when the fields have rolled down to their 
minima, that is at (/> = and a = 62/(264) for a positive initial condition for a. Then /nl 
becomes 

6 „ „ w(l - 6u;al) 



NL 



8- 



:2(l-2^a2)^ 



(5.4) 



where uj = 64/62- Note that within our approximation /nl depends only on the value of cr* 
at horizon crossing once we have fixed the ratio oj. Since the turning of the fields occurs 
only a few e-folds before the end of inflation, we will explicitly assume ~ C/ is a good 
approximation for nearly all the period of inflation. Then we can solve the field equations in 
the slow-roll approximation to find 
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where r = 62/02. 

The time of horizon crossing 



time tfir, 



tfin ~ 60 can be approximately found from the final 



0/4 and thus we calculate the values of the fields at horizon exit as functions 



of the initial conditions and <jq. Using these results in /nl we find 
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(5.6) 



where (/>o = 0o/(2\/6C)/^) 
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We now check the dependence of the above expression on the initial condition ctq. Since 
we assumed that | eio-* | ^ 1 and 1^ ~ ?7 we examine the case ao <^ 1 where 



|/NL = $(l-2u.ag<Ag'^) (5.7) 

rv 



up to second order with respect to do. The parameter /nl becomes maximal if uj 
(^g (4crg) and its value is then 



l-2r 



^ f _ V^O fr 



Since (j)o > 1, the smaller the ratio r and the smaller the initial value of the field a, the higher 
is the value of /nl ■ 

Nevertheless one has to assure that the turn of the field trajectory does not occur too 
late (too close to the end of inflation) , so that the isocurvature mode will have had the time 
to disappear before the end of inflation (so that we can directly extrapolate the results at the 
end of inflation to the time of recombination and do not have to take further evolutionary 
effects into account) and the oscillations of the heavy field do not coincide with those of the 
light field. The higher is the ratio w, the larger is /nl, but then the minimum of the potential 
approaches fio and consequently there is less time available for vi2 and thus for the adiabatic 
perturbation to become constant. This turns out to be a non-trivial requirement: although 
we do not claim to have scanned the whole parameter space of the model, we could not find 
parameter values that passed the above test and at the same time yielded a very large /nl- 
The values we have chosen to work with in this paper respect the above condition and using 
expression (5.7) we expect to find — (6/5)/nl ~ 2. 

If one were to take 64 = 5 instead of 2, one would find — (6/5)/nl ~ 4, but in that case 
the turn of the fields occurs too near the end of inflation so that the isocurvature mode will 
not have disappeared completely by the end of inflation. Looking at the contributions of giso 
and gint separately, we see that even in that case gint has already gone to a constant while 
giso is still decreasing towards zero, so that we feel reasonably confldent that the estimate is 
good even for that model, but we cannot be absolutely certain without a better treatment 
of the end of inflation, which is beyond the scope of the present paper. 

In figure 8 we plot /nl for the model (5.3) with the parameter values described above. 
Again a notable feature comes up during the turn of the fields. It comes from the isocurvature 
term of (4.39) that gets very big during the turn of the field trajectory, but as soon as 
the fields relax it vanishes again. We do not plot gsr separately since it turns out to be 
negligible. Note how the flnal value of /nl depends only on the integrated effect, as the 
isocurvature contribution has vanished. The final slow-roll analytical value is calculated to 
be — (6/5)/NL,sr = 2.15 while the values obtained numerically by our formalism and the 6N 
formalism are — (6/5)/nl = 1-43 and — (6/5)/NL,5Af = 1.48, respectively. We see excellent 
agreement between the exact numerical result of our formalism and the 5N one. The slow- 
roll analytical result does very badly during the turn of the field trajectory, when slow roll 
is badly broken, but gives a reasonable estimate (within 50%) of the final value. 

Finally in figure 9 we plot the spectral index for this model. Its value is in the range 
of the 68% confidence levels after 7 years of WMAP observations [1] , but lies near the upper 
limit. 
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Figure 8. In the first plot we sliow tiie non-Gaussianity parameter /nl as calculated in our formalism 
exactly (red line) and witliin tlie analytical slow-roll approximation (blue dot-dashed line) as well as 
numerically in the context of the 6N formalism (black dashed line) for the model (5.3). In the second 
plot we show again the total /nLj now as the black solid line, and split it up into the isocurvature 
contribution proportional to giso (red dashed line) and the integral contribution proportional to gint 
(blue dot-dashed line). We use the same model as in figure 6. 
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Figure 9. The exact numerical spectral index for the same model as in figure 6. 



6 Conclusions 



The study of the non-Gaussianity produced by inflation models has become a hot topic of re- 
search, since the recent observations of WMAP and in particular the imminent ones of Planck 
will allow us to constrain and discriminate inflation models based on their non-Gaussian 
predictions. In this paper we investigated the super-horizon bispectral non-Gaussianity pro- 
duced by two-field inflation models. To this end we further worked out the long- wavelength 
formalism developed by Rigopoulos, Shellard, and Van Tent (RSvT) [6, 12, 13, 76]. 

We derived an exact result for the bispectrum parameter /nl produced on super-horizon 
scales for any two-field inflation model with canonical kinetic terms, equation (3.6). The re- 
sult is expressed in terms of the linear perturbation solutions and slow-roll parameters. How- 
ever, no slow-roll approximation has been assumed on super-horizon scales, these parameters 
should be viewed as short-hand notation and can be large. In particular this means that the 
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result is valid for models where the field trajectory makes a sharp turn in field space so that 
slow roll is temporarily broken. On the other hand, we did assume slow roll to be valid at 
horizon crossing in order to remove any dependence on the window function and to use the 
analytic solutions for the linear mode functions. Observations of the scalar spectral index 
seem to indicate that this is a good approximation. Note that the assumption of canoni- 
cal kinetic terms is not a fundamental one: the basic equations of the formalism of RSvT 
are given for more general kinetic terms. We just did not want to complicate the notation 
and expressions in this paper with the covariant derivatives and additional curvature terms 
needed to treat the general case. 

The result can be split into the sum of three parts, multiplied by an overall factor (except 
for a small slow-roll suppressed term that is the single-field contribution produced at horizon 
crossing). This overall factor is proportional to the contribution of the isocurvature mode to 
the adiabatic mode, which is only non-zero for a truly multiple-field model where the field 
trajectory makes a turn in field space, as parametrized by a non-zero value of the slow-roll 
parameter rj-^. (Effectively) single-field models do not produce any non-Gaussianity on super- 
horizon scales, since the adiabatic perturbation is conserved in that case. The three parts in 
the sum are: 1) a part that only involves slow-roll parameters evaluated at horizon-crossing 
and hence is always small; 2) a part proportional to the pure isocurvature mode; and 3) 
an integral involving terms proportional to the pure isocurvature mode. Since the adiabatic 
mode is not necessarily constant in the presence of isocurvature modes, we only consider 
models where the isocurvature mode has disappeared by the end of inflation, so that we 
can directly extrapolate our result at the end of inflation to recombination and observations 
of the CMB. However, this automatically means that the part 2), although varying wildly 
during the turn of the field trajectory, cannot give any persistent non-Gaussianity that can 
be observed in the CMB. This means that any large non-Gaussianity on super- horizon scales 
in models satisfying this condition will have to come from the integrated effect in part 3). 

The exact equation (3.6) is the basis of our numerical studies. However, to gain further 
insight we tried to work out the integral analytically. For this it turns out that the slow-roll 
approximation is necessary. Even then the integral can only be done explicitly for certain 
specific classes of potentials, among which are product potentials, W{(p,a) = U{(l))V{a), and 
generalized sum potentials, W{(j),a) = {U{(p) + V{a))'^ . We found that, with our assumptions 
on the disappearance of the isocurvature mode, no product potential can give large non- 
Gaussianity, nor can any simple sum potential with equal powers, W{(l),a) = a4>^ + /3cj^. 
However, we found conditions under which the (generalized) sum potential can give large 
non-Gaussianity (here defined as /nl larger than unity), and we have described an explicit, 
simple model that does. It consists of a heavy field rolling down a quadratic potential while 
a light field sits near the local maximum of a double-well potential. When the heavy field 
reaches zero and starts oscillating, the light field takes over and rolls down, so that there is a 
turn of the field trajectory in field space. We studied this model numerically, using the exact 
results, to confirm our analytical predictions. 

In deriving equation (3.6) we assumed that all three scales cross the horizon at the same 
moment. However, this is not a necessary assumption, and we also generalized the result 
to an arbitrary momentum configuration. We find that going to the squeezed limit, where 
one of the momenta is much smaller than the other two, even when remaining within the 
resolution of the Planck satellite {k' ~ lOOO/c), the result for /nl can be reduced by about 
10%, depending on the model. We stress that we are discussing /nl here, so this effect is 
unrelated to the well-known result that the local bispectrum peaks on squeezed momentum 
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configurations, which is due to the momentum behaviour of the power spectrum, which has 
been divided out in /nl- However, exactly because of this latter effect, the squeezed limit is 
very relevant for the computation of /nl- 

We have worked out and included the second-order source term at horizon crossing in 
the long-wavelength formalism of RSvT, a contribution that had been missing so far. This 
is the only change of the basic formalism with respect to the paper [6] by RSvT. While this 
additional term is always small for the models we consider and hence numerically insignificant, 
from an analytical point of view it means we could now compare our results directly to the 
so-called as defined in the 6N formalism [19, 62, 63]. Some of the potentials we studied 
had already been worked out using that formalism and where available we compared our 
analytical results and found perfect agreement. We also compared our exact numerical results 
with those obtained using a numerical 6N treatment for models where slow roll is broken 
and the analytic results cannot be trusted, and again we found excellent agreement. 

We showed that the long-wavelength formalism of RSvT represents a viable alterna- 
tive to the 6N formalism to compute the super-horizon non-Gaussianity produced during 
inflation, allowing us to obtain and verify results in a different way. Moreover, the long- 
wavelength formalism has a number of advantages that can make it preferable in certain 
situations. Very importantly, our formalism allows for a simple physical interpretation of 
the different parts in terms of adiabatic and isocurvature modes, providing insight into the 
behaviour of the different transient and persistent contributions to /nl- While we did not 
pursue this in the present paper, the formalism also provides the solution for the second- 
order isocurvature perturbation and hence the isocurvature bispectrum could be computed 
as easily as the adiabatic one. 

From our studies it has become clear that the condition on the disappearance of the 
isocurvature mode by the end of inflation is a very strong constraint. It significantly reduces 
the possibilities for a large, observable value of /nl produced during inflation. Note, however, 
that we chose to impose this condition only to be able to neglect the further evolution of 
the adiabatic mode after inflation; it is in no way a necessary condition for our formalism 
during inflation. In future work we would like to relax this condition, which means that the 
adiabatic mode would no longer necessarily be constant after inflation, and hence will require 
a much better description and understanding of the evolution of the perturbations during the 
transition at the end of inflation and the subsequent period of (p)reheating. In conclusion, 
while a lot of progress has been made over the past few years regarding the non-Gaussianity 
produced in multiple-field inflation, more work still remains to be done. 
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A Basis improvements 

In [8, 77] an orthonormal basis in field space was introduced, with substantial advantages for 
calculating and interpreting quantities in multiple-field inflation (see also [78] for a special 
two-field case of this basis). The basis was defined as follows (note that e.g. 11 is the vector 
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containing components 11^). The first basis vector ei is the unit vector in the direction of 
the field velocity. Next, the direction of the basis vector 62 is given by the direction of that 
part of the field acceleration that is perpendicular to ei. This orthogonalization process is 
then continued with higher-order time derivatives, until a complete basis is found. Defining 
the generalized rj parameter as 

(n) ^ ^dtr-'u 

' //"-in ^ ^ 

(with 77 = 77^^) and ^ = ^7^^^), we can then define the basis vectors via an iterative procedure 
as 

(n) _y-n-l (n) 

rin 

^ ^ = p. .r,(")_ 

basis vectors, 

to be non-negative: 



for n > 2, with ei = Il/n and ry) ' = ei ■ 77'-"-'. Basically there is an arbitrariness in the 

(n) 

choice of sign of the basis vectors, which in the original definition was fixed by choosing r/n 



n-l 



(old definition) (A.3) 

=1 ' 

While being a perfectly valid choice analytically, this choice does mean that certain basis 
vector components and slow-roll parameters make sudden sign flips when one or more fields 
are oscillating, and that is hard to deal with numerically. Hence we now propose a different 
choice for 77^"^ , which is identical except for the overall sign, and which eliminates the sudden 
sign flips: 

?7^") = -eAi...A„efi • • • e^Ti'v^''''^", (new definition) (A.4) 

where e is the fully antisymmetric symbol. Prom the fact that rj = Y17=i ''li^^^i it immediately 
follows that 

£Al...A„e^•••e^ = -l, (A.5) 

so that this choice means that the basis has a definite handedness. Note that in the case 
where the fields do not oscillate, the two definitions have the same overall sign (hence the 
choice of the minus sign) . To have the expressions for the time derivative of the basis vectors 

(n) 

and the rjn unchanged, we see that we also need the relation 

to be satisfied. Then all results and expressions developed with this basis are unchanged. 
An interesting consequence of these relations, including the orthogonality relation 

is that for the cases of two and of three fields we have sufficient conditions to write all basis 
vectors in terms of ei, without knowing anything about the dynamics. For two fields we have 

62 = {el -e{), (A.8) 

with (e})^ -h (ef)^ = 1, and for three fields 

62 = {el + el -e\ + el -e\ - e?), 

ez = {\- {e\f - e\el -\ + {e\f - e\el \ - {eff - e}e?) , (A.9) 
with {elf + (ef)2 -|- (ef)^ = 1 and e^ef - e}ef + ejef = -5. 
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B Computation of the second-order source term 



(2'] 

To compute the second-order source term h)2' of (2.13) we choose a gauge characterized by 
the constraint ei^c/?"^ = 0, so that d = a and C2 = —{f^/^/^)e2A'^^, where we have spht 
up (j)"^{t^x) = 4>'^{t) + ip^{t,x) and lna{t,x) = Ina(t) + a{t,x). Note that we use here 
subscripts 1 and 2 to indicate the adiabatic and isocurvature components of C (instead of 
superscripts). On long wavelengths this gauge reduces to the uniform energy density gauge. 
We will drop the tilde (that we introduced in section 2.3 to denote quantities in the uniform 
energy density gauge) and omit the superscript on first-order quantities in this appendix 
in order to lighten the notation. 

We can use the momentum constraint (the equivalent of (2.3) in the proper gauge) to 
solve for the lapse function to first order by setting = 1 -|- A^i, and we find A'^i = C/H. 
Following Maldacena [5] we first write the quadratic action in the ADM formalism. For that 
we only need to compute A^i to first order. We find that to leading order in slow roll the 
action takes the form 

S2 = I dtSxL2 = I dtd^x a^ejcf + C| - Vi/CiC2 + 2x^^6(2 1 ■ (B.l) 

To obtain this expression we used the gauge constraint to simplify the integrand as well as 
(2.4) and (2.5). 

For the third-order action we find for the leading-order slow-roll terms [69] 
83 = J dtd^x I a3e2 (interaction terms) - ^ + (e + r/ll)CiC2 + ^ + ^C?^ 

1^-^ - T + - ^ ^^^^ + ^ j I' ^^-2) 

where the auxiliary quantities Qi are defined as 

Qi = --^eiA^f^y (B.3) 

and the last two terms of the action are proportional to the second-order equations of motion. 
In order for them to vanish we perform a redefinition of Ci and 

Ci = Cic - ^ - + Ci - ^ C1C2 + 

C2=C2e-f + (. + .ll)ClC2 + ^ + ^C?. (B.4) 

The terms Qi can be rewritten using the definitions of the gradients of the curvature pertur- 
bations at second order. While working in the uniform energy density gauge we can use the 
equivalent of (2.3) to find the following constraints to first and second order [69]: 



H(j)Adnp(^) = 0, 



'Qi = ^Cl + #-5-5^fe.C2), (B.5) 
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while expanding (^2i up to second order in the same gauge gives 



1^(2) 



a,; 



+ ^C2Cl 



Then the redefinitions of the gradients become 



Ci + 



(2) 



Clc + 



CiCi _ C2C: 

H H 



' + ^ (C? - Cl) - r/^CiC2 + 5-^9^ ( 



(B.6) 



(B.7) 



C2 + 



(2) 



> , C2C1 , C1C2 
C2c+^ + ^ 



+ V (^1 - ^2) + + ^")CiC2 - 5-29^ ( |5.Ci ) . (B.8) 



Inspecting (2.30) and (B.6) we see that in the uniform energy density gauge (mi = diQm both 
for the adiabatic and the isocurvature component. Since we perform our main calculation in 
the flat gauge and we use the variable Qmi rather than C^, we want to transform the above 
redefinitions to this gauge by the simple gauge transformation (C.3): 



(2) _ ?(2) Cl 



H 



(B.9) 



where we used again the tilde to denote the uniform energy density gauge. We find 



a(2) 



di 



Clc + 



CiCi _ C2C: 

H H 



2 ^ e + 7? 



Cl 



did 



C2i + 



A2) 
^2i 



' , C2C1 , C1C2 , r]^ 

C2c + -^ + ^ + Y 



(Cf-cl) + (e + r?")CiC2-a-2a^(|9X: 



C2 



(B.IO) 

Cl O A 

- ^5.C2 
(B.ll) 

Note that after horizon exit Ci = 2Hr]-^(2 and C2 = —HxC2 (the latter is valid under the 
slow-roll approximation only) in the gauge used in this appendix, so that the expressions 
simplify. For a field redefinition of the form C = Cc + ^Cc (note that in the equations above 
we have not added the subscript c explicitly in the quadratic terms, since to second order it 
makes no difference) the three-point correlation function can be written as 



(CCC) = (CcCcCc) + 2A[(CcCc)(CcCc) + cyclic] 



(B.12) 



Hence the elements of Liab and Ni^b are just the coefficients of the various products of Ci 
and C2 in the redefinition of Ci multiplied by 2, and similarly for L2ab and A^2a6- Note that 
the local terms Labc correspond to the terms between the square brackets, and the non-local 
terms Nabc to the terms outside. This leads to the explicit slow-roll expressions given at the 
end of section 2.2. 



C Gauge issues 

C.l Gauge invariant quantities 

The use of spatial gradients was first advocated in [79] in the context of the covariant for- 
malism. Later on the authors of [12] constructed invariant quantities under long-wavelength 
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changes of time-slicing, by considering the following combination of the spatial gradients of 
two spacetime scalars A and B: 



Ci = diA-^diB, (C.l) 
and formed among others the quantity 

Ci = d,a - ^dip, (C.2) 

OtP 

which when linearized is just the gradient of the comoving curvature perturbation ^. Here 
a is the logarithm of the scale factor and p the energy density. In [55] it was argued that 
this particular combination could give rise to a second-order gauge-invariant quantity ^^^^ as 
follows: 

= Cf ^ - ^Cf \ (C.3) 

where C*-^^ can be related to yet another gauge-invariant second-order quantity defined by 
Malik and Wands in [53], 

^(2)^^(2) _^(1)2 

(2)1 

We can easily show that (C.3) is equivalent to taking the gauge transformation of Q 
from homogeneous expansion time slices to uniform energy density time slices. Denoting as 
tilded the quantities in the new time slices, where T = t + At (no space transformation is 
required), we can write 

r(l)l ^ (1)1 



Cj^)^ = C^^)^ - At (C.5) 



The quantity At is the time difference between the two time slices and can be evaluated by 
either comparing a (see also (2.31)) or p in the two gauges: 

p{t, x) = p{t + At) = p{t) + At ^(t), (C.6) 
which expanded to first order becomes 

pW(t)+p«(t,x)=p(t) + At^(t). (C.7) 
Hence we find that p^'^^ = p and 

At = (C.8) 

~(2)l 

We conclude that our definition of the second order Q in the uniform energy density gauge 
agrees with the gauge- invariant quantities defined by other authors. 
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C. 2 Gradients and locality 

As a consistency check we want to verify that Q is indeed a total gradient, as it should be 
according to the first part of this appendix and (2.30). Taking expression (2.19) (corrected 
by the gauge transformation), 

Cl^^^ = -{diVe,)vf, |- 2r]^Gif{t, U)G2e{t, U) + dt'Gia{t, t')AabcGbe{t', U)G^f{t\ t,)| 

+Gia{t,U)Laef*di{Ve*Vf^) + Gia{t,U)Naef*{diVe*)vf^, (C.9) 

and rewriting it using (D.l) we find 

Qf'^^ = -{diVe.)Vf, ^ab{U)Gia{t, U)6beSlf + dt'Gia{t, t')AabcGbe{t\ U)Gsf{t' , U) 

— J^dt' AabAicGia{t,t')Gbe{t' ,t^)Gcf{t' ,t^,) — Gia(t, t*) A^ae/* ^ 
+Gla{t,U)Laef*di{Ve*Vf^). (C.IO) 

This expression should be symmetrical under the interchange of the indices e and /. Notice 
that the last term is automatically symmetrical. 

The anti-symmetrical part of the two integrands turns out to be proportional to 

Ta = G'23(t', t*)G32(t', t*) — G22(i', i*)G33(t', t*). (CH) 

We explicitly check the exact numerical value of this quantity and find it to be zero. For 
this it is crucial that we have defined as the time a few (about 3) e-folds after horizon 
crossing. The reason is that the long-wavelength approximation we use in all our derivations 
is only valid once the rapidly decaying mode can be neglected, which takes a few e-folds. If 
the above quantity were to be evaluated before that time, it would not yet be zero. Note 
that in the slow-roll case Ta is identically zero according to (4.2), since within the slow-roll 
approximation the decaying mode is neglected by construction. The above means that Cl 
is well defined only after it is well outside the horizon, where we can neglect the decaying 
mode. The case where the decaying mode can remain important is treated in the one-field 
case in [68]. 

The remaining non-integral terms between the braces can be explicitly checked to cancel 
when taking the slow-roll limit at horizon crossing: the first term of the first line of (C.IO) 
gives 

- {diVe*)vf^Aab{t^)Gia{t,t^)5beSlf = (2r]^ - xGl2{t,U)^ ivudiV2* - V2*diVu) , (C.12) 

while the terms arising from the non-local contribution N^ef* gives exactly the same but with 
opposite sign. Hence we see that within the conditions of the long- wavelength approximation 

~(2)1 

Q is indeed a total gradient, as it should be. 
D Detailed calculations 

D. l Relation between space and time derivatives 

We begin by proving that 

Aabl = -^dtAab (D.l) 
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for all gauges with dilna = 0. Actually the statement is more general: the prefactor of 
in the expression for dif, where / is any function of H,cf),Tl, is equal to —{dtf)/{NH) for 
gauges satisfying dilna = 0. 

We start by computing the time derivative of /: 

= ^ {dnf m + V,/ . 9,0 + Vn/ • d,U) 

= -HedHf + —ei-V^f + ^Hrj-Vnf, (D.2) 

K K 

where we used the definitions of the slow-roll parameters e and ij, and of the basis vector ei. 
On the other hand, the spatial derivative of /, in a c?j Ina = gauge, is given by 

dif{H, ct>, U) = dnf d,H + V^/ • 0,0 + Vn/ • diU 



HedHfei-C,-—Ci-^4>f (D.3) 

Kj 



2e 



H{^e, + {e + - e(ei • Ci)ei ■ Vn/, 



NH 

using the constraint relations for diH, di(f), and 9,11 given in [6, 12]. (Note that some time 
and space derivatives have to be replaced with their covariant (in field space) version to 
make contact with the more general expressions given in those papers that take into account 
a non-trivial field metric.) Taking the components of ^ and 6 in the basis defined in A, not 
forgetting the relation • 0i = + N H ZmnCi , with Z21 = rj-^ [6, 12], we prove the stated 
relation, of which (D.l) is a special case. 

D.2 Derivation of equation (3.6) 

In this appendix we work out the last term of (3.2), which has to be added to the result for 
the second term (/) given in (3.5), to derive the final expression (3.6) for /nl- We call the 
sum of these two terms J: 



J _ I 



+ / dt'GiaAabcVbmVcn, (DA) 



which is 

-'5m2<5n2 2r?;;fx* Gi2.{t,U) (D.5) 
+<5™2<5„2 /,* dt' [ (^122 - 4(r?^)2 + A322G13) {V22? + (^33 + 2r/^) Gi^{v^2? 



+ (^123 - 4r/^Gi2 + (^323 + ^32 + 2r/^^33 + 2??^)Gl3) ^^22^^32 



We remind the reader that the bar on top of an index (c) means that it does not take the 
value 1 and that a subscript * means that a quantity is evaluated at t^. The explicit form of 
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the matrix A is given in [6]: 

ii2i = 2e7?^ -47?ll7?^ + 2^^, 

Ai22 = -6x-2e7?ll -2(7?ll)2-2(?7^)2, 

Ai2z = -6-2r?", 

^321 = -12e7/ll - 12(r/-^)2 - 6ex - 8e^ - 2{)e^r{^\ - Aeirj^f - I2e{r]^ f 

+16r/ll (77^)2 - ee^ll - I2r]^i^ + 3(VKiii - 1^221), 
^322 = -24e7/-^ - 12r/llr/-^ + 24r/-^x - 12e%-^ + 8(r/ll)2r/^ + 8(77^)^ 

-8ee^ - 477!'^^ + 3(1^211 - ^-222), 
^23 = 12r?^ - 4er/^ + Sri^r]^ - 4^^, 
^331 = -262-4er/ll+2(r?ll)2- 2(77^)2 -2^11, 
^332 = -4677^ -2^^, 

^333 = -2r/^, (D.6) 
while the rest of the matrix elements are zero. Using these expressions we have 

^33 + 2r/^ = 0, 

^23 + ^32 + 2r/^^33 + 2??^ = 1877^ - 47?^, 

ii23 = -2^33 + 2e + 2r/ll, 
^22 - 4(7?^)2 = -2A32 + 2e + 27)11 , (D.7) 



(5m2^ni - 2ri^+ x*Gi2(t, U) + ^32*Gi3(t, U) - x*Mz*Giz{t, hfj 



so that we can write 
J 

-<^m2'5„2 2ri^x* Gisit, t*) 

+Sm2Sn2 fl dt'[ 2fe)'^(e + V^^) + 2^22^1^32 " 4 (77^^12 + 77^^13) ^^(^22)' 

+2(e + r/l')-i;22W32 + ^322Gl3(i;22)^ + 187/-^Gl3t;22l^32 • (D.8) 

Doing integrations by parts on the three terms in the third line we obtain 

^ = Sm2Kl (-2r]^+ X*Gl2it, U) + A32*Gl3(t, U) - X*^33*G'i3(t, U)\ 

+2(5.m25n2^ -rj* X*Gi3{t, U) - (e* +77!!) + x* + 'n^Gi2{t,U) 
+V*Gi3{t, Q + (e + 77ll)(i;22)^ + ^^22^)32^ 

+26^m2Sn2 J^dt' - 2{ri^)\v22f - {e + r]^^)v22V32 - {V32f + 9V^Gi3V22V32 

+1 [A322 + 2ii^ + 277^^33 + 277^^32) Gi3(7;22)^] . (D.9) 

The following relation (derived by taking two time derivatives of the field equation) can 
be used to remove higher-order slow-roll parameters: 

(4) / j| \ 

W^u = - 1 - |- Urn + (2e + 7yll)7?™ + eri^mi - n^W^2. (D.IO) 
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Explicitly, for m = 1 and m = 2 in the case of two fields, this becomes 



Will = -^r? 



(4)1 



1^211 = -^r/ 



(4)- 



(D.ll) 



Using the second of these relations, as well as the explicit expression for ^322, we find that 

^322 + 2f?^ + 27)^^33 + 2r/^^32 (D.12) 
= 2Ar,^x - 12??llr/^ + Ui'q^^r]^ + I2{ri^f - V^" - 47?"^^ - 3(#2ii + ^^222)- 

We now drop boundary terms that are second order in the slow-roll parameters at horizon 
crossing, since it would be inconsistent to include them given that the linear solutions used 
at horizon crossing are only given up to first order. Then the result is 



J 



Sm2Sni{-2r]^ + X* V12) + 2dm2Sn2 (^-e^ - r]l + X* + vivi2 + (c + r?")(w22)^ + f 221^32) 



+25m2'5n2 //^ dt' 



-2(?7-^)2(t!22)^ - (e + T)'"22V?,2 - {Vmf + 9rj-^Gi3V22V32 



--(VF2II+W222) Gl3(^22)^ 



(D.13) 



Inserting this into (3.1) gives the final result for /nl in (3.6). 
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